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Abstract. We examine the time evolution of cold atoms (impurities) interacting with an 
environment consisting of a degenerate bosonic quantum gas. The impurity atoms differ from the 
environment atoms, being of a different species. This allows one to superimpose two independent 
trapping potentials, each being effective only on one atomic kind, while transparent to the other. 
When the environment is homogeneous and the impurities are confined in a potential consisting 
of a set of double wells, the system can be described in terms of an effective spin-boson model, 
where the occupation of the left or right well of each site represents the two (pseudo)-spin states. 
The irreversible dynamics of such system is here studied exactly, i.e., not in terms of a Markovian 
master equation. The dynamics of one and two impurities is remarkably different in respect of the 
standard decoherence of the spin - boson system. In particular we show: (i) the appearance of 
coherence oscillations, (it) the presence of super and sub decoherent states which differ from the 
standard ones of the spin boson model, and (m) the persistence of coherence in the system at long 
times. We show that this behaviour is due to the fact that the pseudospins have an internal spatial 
structure. We argue that collective decoherence also prompts information about the correlation 
length of the environment. In a one dimensional configuration one can change even stronger the 
qualitative behaviour of the dephasing just by tuning the interaction of the bath. 



Collective decoherence of cold atoms coupled to a Bose-Einstein condensate 



2 



1. Introduction 

The reasons of the great interest for the physics of ultracold atoms in recent years are manifold. On 
the one hand experimentalists have reached an unprecedented control over the many-body atomic 
state with very stable optical potentials and by the use of Feshbach resonances which allow one 
to change the scattering length of the atoms [lj. In this context the tremendous experimental 
results that have been achieved include: the observation of the superfluid-Mott insulator transition 
for bosons [2], one dimensional strongly interacting bosons in the Tonks-Girardeau regime [3] and 
Anderson localization [H [5] . On the other hand new experimental challenges come from different 
theoretical proposals for using this system for quantum information processing [6] and as a quantum 
simulator of condensed matter models (see for example [3, El [9] and references therein). 

Not only can ultracold atoms simulate Hamiltonian systems, but such systems also offer a way 
to engineer non classical environments. Thanks to the flexibility of quantum gases, a broad range of 
regimes of irreversible dynamics of open quantum systems and in particular of spin-boson systems 
can be explored [TOj El H21 H21 M ■ 

In the present paper we propose a new way in which an instance of the spin - boson model [15] 
can be realized with a suitable arrangement of interacting cold atoms. In particular we analyse a 
system consisting of cold impurity atoms interacting with a degenerate quantum gas of a different 
atomic species. This setup makes possible the superposition of two independent trapping potentials, 
each being effective on one atomic species only, while transparent to the other. When the quantum 
gas is homogeneous and the impurities are confined in a potential composed of double wells, the 
system can be described in terms of an effective spin-boson model, where occupations of the left 
or right well represent the two (pseudo)-spin states. At variance with other setups, where the role 
of the pseudospin is played by the presence or absence of one particle in a trapping well [16], by 
the vibrational modes of a single well [17] or by internal electronic levels [H], in our case each 



pseudospin has a spatial dimension, namely the separation between the two minima of the impurity 
double well. This introduces an effective suppression of the decoherence due to low frequency modes 
of the environment and leads to unusual and interesting phenomena, like oscillations of coherence 
at finite times and the survival of coherence at long times. Further novel features appear when one 
considers the irreversible collective decoherence of a systems of two impurities. In this case we still 
predict the existence of subdecoherent and of superdecoherent state, but with the interesting fact 
that their role is exactly the opposite from what one observes in conventional spin-boson systems. 
Further interesting features appear when one considers how the collective decoherence rates change 
as a function of the impurities' separation and the effects of dimensionality of the system. 

In discussing our investigations, for the sake of simplicity we shall consider an experimental 
setup where the impurity atoms are trapped by a periodic (optical) lattice. We like to stress, 
however, that our findings do not depend on the lattice properties (e.g., periodicity) but for the 
numerical results. Other setups, such as microtraps on atom chips or quantum dots, just to mention 
a few, can be equally envisaged. 
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Figure 1. A Bose-Einstein condensate (yellow region) confined in a shallow harmonic trap Vb(x) 
interacts with cold impurity atoms each of which is trapped in a double well potential Va(x) (gray 
circle). The distance between two wells in the same trap is 2L and the distance between adjacent 
traps is 2D. 

2. The Hamiltonian 

Our system is composed of a cold quantum gas of bosonic atoms and a sample of cold atoms 
separated from each other and immersed in the quantum gas. In presenting our investigations, we 
shall use the words 'reservoir', 'bath' and 'environment' as synonyms to indicate the quantum gas, 
since its properties are not the focus of the present paper. 

The second-quantized form of the Hamiltonian of the impurities+bath system takes the form 
(see also Ref.pSJ) 

H = H A + H B + H AB (1) 

where 



x 



2mA 



¥(x) 



(2) 



is the Hamiltonian of atomic impurities, described by the field operator ^(x) in the trapping 
potential Va(x) which creates a set of double wells of size 2L and separated by a distance 2D, see 
Fig.IH 



'd 3 x $t(3 



+ VB(x) + ^$t( x )$( x ) 



2m b 



$(x) 



(3) 



is the Hamiltonian of the bath, composed of N ^> 1 bosons, represented by the field operator $(x) 
and confined by a trapping potential V#(x) and gs = 47rft 2 ag/m^ is the boson-boson coupling 
constant, with as the scattering length of the condensate atoms, and 
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#ab = 9ab Jd 3 x *t(x)&+(x)d(x)*(x) (4) 

describes the interactions between the impurities and the bath; here Qab = ^T^O'Ab/^ab is the 
coupling constant of impurities-gas interaction, with a^s the scattering length of the impurities-gas 
collisions and ttiab = fnA m B/( m A + itib) their reduced mass. Both impurity and bath atoms are 
described in the second-quantized formalism. The field operator of the atomic impurities 



£(x) = £a i)p¥ > i)P (x) (5) 

can be decomposed in terms of the real eigenstates (fi )P (pc) of impurity atoms localized on the double 
well % of the potential Va{x) in the p th state, with energy Tiuo i p and the corresponding annihilation 
operator a i . We assume that the wavefunctions of different double wells have a negligible common 
support, i.e., y? f (x)^?,y ini (x) ~ at any position x. 

We treat the gas of bosons following Bogoliubov's approach (see, for instance, Ref.[T9j) and 
assuming a very shallow trapping potential V B (x), such that the bosonic gas can be considered 
homogeneous. In the degenerate regime the bosonic field can be decomposed as 



$(x) = 0Vo $ (x) + <f$(x) = y/iVo $o(x) + £ («k(x)c k - <(x)4) (6) 

k 

where $o( x ) is the condensate wave function (or order parameter), N < N is the number of atoms 
in the condensate and the annihilation and creation operators of the Bogoliubov modes 

with momentum k. For a homogeneous condensate $o( x ) = V being the volume. Its 

Bogoliubov modes 



1 1 (tk + n g B T e ik " 

= V 2 1" & — V W (8) 



have energy 



2 -| 1/2 

2e k n g B + e k , (9) 



where e k = ^V/(2mg) and n = N /V is the condensate density. As one can see from (jUj), low- 
energy excitations have phonon-like (wave-like) spectrum, whereas high-energy excitations have 
particle-like spectrum. The condition for wave-like excitations is e k <C n gB, i.e., k <C Ay/WnoOs, or 



equivalently <C 2mBC s /h, where c s = JnogB/rn is the speed of sound at zero temperature. Note 
that |w k | = 1/yV and |t> k | = describe the limiting case of iV 1 non-interacting bosons, each 
with energy i? k = e k . 
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Inserting Eqs.fjS]) and (jSJ) into the Hamiltonian flTJ we get 



Ha = ^2 h^hp 



^i,p^i,p 



for the impurities, 



(10) 



Hb — Hcond + -f^Bog 

for the quantum gas, with 



(11) 



#Cond = N Jd*X $*(x) 



for the condensate and 



*o(x) 



(12) 



Bog 



for the collective excitations (Bogoliubov modes) of energy E k in the condensate, and 
H A b = 9ab °\,p% q N o / rf3 ^, P ( x )^i,g( x )l $ o(x)| 2 + 

i PA V J 

+ \/^0^ a k /^ 3 ^,p(x)^, g (x) ($*(x)« k (x) - $ (x)w k (x)) 
k J 

+ / d 3 x^, p (x)<^(x) ($ (x)<(x) - $*(xK(x)) 



(13) 



(14) 



for the interaction Hamiltonian; the terms which are quadratic in the Bogoliubov excitation 
operators c, give negligible contributions and have been omitted. The first term in (|14p describes 
transitions between impurities' vibrational states due to the condensate, whereas the remaining 
terms describe similar transitions induced by the collective excitations in the condensate. In a 
homogeneous condensate, transitions between different vibrational eigenstates of the impurities 
induced by the condensate are suppressed, while all vibrational states <^i, p (x) get an energy shift 
Sid, 



g No d 6 x |$ | 2 (x)^ iP (x)^(x) 





AH c 

n g = duji t p 



for p ^ q 
for p = q 



(15) 



so the contribution of the first term in (TBI) can be included in the definition of u^p. 

In the limit of deep, symmetric wells in each double well and separated by a high energy barrier, 
the tunneling between adjacent wells is suppressed. In this regime the ground states (p^z and <p i; R 
of, respectively, the left and right well of double well i are well separated in space with vanishing 
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spatial overlap, their coupling to the excited states becomes negligible and the total Hamiltonian 
further simplifies into 



H 



E E *»t*K + E eAK + E a E E K>A + 



(16) 



i, p=L,R k 

where we have defined the coupling frequencies 



p=L,R k 



tf*k = (l«k| - kl) / d 3 * K P (x)|V k ' x (17) 

and rip = a\ p a i p is the number operator of impurities in the double well i in the well p = L,R. 

We consider the case where each double well is occupied by at most one impurity atom. This 
allows us to describe the occupation of the left and right well of each site in terms of pseudospin 
states. Introducing the Pauli operators as h l L — (1 — a*)/2, h l R — (1 + o~l)/2, the Hamiltonian ([TBI) 
takes the form of the independent boson model [20] 



H 




E 



L,k 



K + E + n 



L,k 



+ 



E Kk - «£k) ^ + E Kk + 



where a constant energy shift has been omitted. We note that spin-boson systems with larger spin 
values can be realized in the same way with higher occupation of the double wells. 

The effects due to quantum noise on coherent superpositions of states of a double well spin- 
boson hamiltonian have been analyzed in the markovian regime. In [2T[ [22[ [23] the effects of a cold 
atom reservoir has been analyzed, while [23] has considered the effects of scattered photons, taking 
into account also the role of the inter-well separation. As we will show in the following section, for 
our system it is possible to carry a full analysis of the impurity dynamics, going beyond the Markov 
approximation. 



3. Exact reduced impurities dynamics 



The dynamics due to the spin-boson Hamiltonian ([TBI) is amenable of an exact analytical solution 
and is characterized by decoherence without dissipation J25[ [26j [27]. The time-evolution operator 



U(t) = exp —iHt/h corresponding to the Hamiltonian ([TBI) can be factorized into a product of 
simpler exponential operators, 



U(t) = exp 



E ^kcjA* 



x exp 



E E Ai(t)&i + a k (t) 4- E E A t(t)^ + <(t) % 



l k 
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x exp 



x exp 



x exp 



ft 2 



4El 



2^ 



ift 2 E/k(*) 



E« (nk,k + ^k) Ej (ng k + » 

4E£ 



L,k 



(19) 



where the functions 
/k(«) : 

a k (t) = 



•E'k, ^k, 

— t — sin — r, 
n n 

h(l- e iE ^ h ) 

2E~k 
h(l- e iE ^ h 



h'.k --l.k j 5 



2E U 



E Kk + n 



2* 



(20) 
(21) 
(22) 



have been introduced for ease of notation. Details of the derivation of f)19p for the time evolution 
operator are given in Appendix A. As in this paper we are interested in the irreversible collective 
decoherence of the impurities we will focus our attention on the conditional displacement operator 



Uo{t) = UUuAt), 



U k , D (,t) = exp 



(23) 



(24) 



Indeed this operator is the one responsible of the decoherence of impurities as it induces 
entanglement between them and the reservoir. Labeling the state of the impurities as = 
\{ni, n 2 , n 3 , . . .}) with n p = 0, 1 denoting the presence of the atom, respectively, in the left or right 
well, the matrix elements of reduced density operator of the impurities are 



P{n p },{m p }(t) = exp [-r { „, i}> .r m j(t)J P{„. p },{m p } (0) 

x exp |z0 

{n p },{m p } 

(t)}exp {t (t)}exp{iA 

{n p },{m p } 

(*)} (25) 

Assuming that each mode of the bosonic environment is in a mixed state pk at equilibrium at 
temperature T the decay exponent contains all the information concerning the time dependence of 
the decoherence process and takes the form 



r {ni},{mj}(*) = ^ 2 E 



E K - ^] (^kk - 



coth 



(26) 



with (3 = 1/K B T. The phase factors <S>{n p },{m p }(t), E {nph{mp} (t) and A {np}){nip} (t), whose specific 
form is given in appendix B, do not play any role in the decoherence [28]. They contain however 
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interesting information on the effective coupling between the pseudospins induced by the consensate 
and will be analysed in a future paper [29J . 



4. Results for the decoherence 



As mentioned in the Introduction, we shall assume that the impurity atoms are trapped by an optical 
(super)lattice, whose form can be controlled and varied in time with great accuracy [301 GHJ- The 
coupling frequencies VL l k are accordingly evaluated in Appendix C assuming an optical lattice, with 
identical, double wells in each site, and deep trapping of impurity atoms in their wells, with identical 
confinement in each direction. Atomic wave functions can then be approximated by harmonic 
oscillator ground states of variance parameter a = yjh/ (mu) [32], where lj is the corresponding 
harmonic frequency. As will be clear shortly, a acts as a natural cutoff parameter, quenching the 
coupling with high frequency modes. 

Specifically, we consider 23 Na impurity atoms trapped in a far-detuned optical lattice and a 87 Rb 
condensate. The condensate density is no = 10 20 m -3 , the lattice wavelength is A = 600 nm, and we 
have taken 2L = A/2 and D = 2L. The depth of the optical lattice is described by the parameter 
a = Vq/Er, Vo being the optical lattice potential maximum intensity and Er = h 2 k 2 /(2m) the 
recoil energy of impurity atoms in the lattice; in our evaluations we put a = 20. Finally, we assume 
a AB = 55a [33], where a is the Bohr radius, for the scattering length of impurities-condensate 
mixtures. This parameter can be modified in laboratory with the help of Feshbach resonances. 



4-1. Single impurity decoherence 

We first examine the decoherence exponent of a single impurity 



r (t) = r {0} , {1} (t)^ 2 £ 

k 



CQS f*) CQth ^k 



El 



n 1 - n 1 

ic R,k "L,k 



(27) 



Such quantity, that will be a useful benchmark in our analysis of the collective decoherence of 
impurity pairs, shows already interesting features. Assuming, from now on, that the condensate is 
at temperature T = 0, we obtain 

r (t) = 8g 2 AB n J2 (Kl - kl) 2 e- fe2CT2 / 2 ^# sin 2 (k • L) (28) 



El 



We note the dependence of T (t) on the length L, where 2L is the distance between two wells 
within each site. The presence of the factor sin 2 (k • L) supresses the effect of the reservoir modes 
at small k. This is clearly understandable: environment modes whose wavelength is longer than 
L cannot "resolve" the spatially separated wells within each site. The consequences of this fact 
will be clear shortly. Replacing the sum over discrete modes to a continuum with the usual rule 
V^ 1 J2k ~^ (2?r)~ 3 / dk, choosing x as azimuthal axis and using well known relations for Bogoliubov 
modes [34], we finally obtain 



rS(*) 



2g 2 AB n o 



oo 
dk 



sin 



2 ^k- 

2h ' 



E k (e k + 2g B n ) 



1 - 



sin 



2kL\ 



2kL J 



(29) 
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Figure 2. To(t) vs. time for a single impurity atom interacting with free bosons (solid line) and 
with a bosonic condensate (dashed line) in three dimensions. The inset shows r (£) for very short 
times < t < 2 /is. 



The superscript c is to remind us that we are dealing with impurities interacting with a condensate. 
For the special case of a bath of noninteracting bosons T^^t) is obtained from f[2"§l) simply imposing 
gs = and E^ = e^. Let us point out that the spectral density, which reads 



has a non trivial form, which at small frequencies, scales as u d+2 for the interacting case, where 
d is the dimensionality of the condensate, and as uj d l 2 for the non interacting case. It is worth 
noticing that while the former case is always superohmic, the latter one is subohmic, ohmic and 
superohmic depending on the dimensionality of the environment. Note that the high power in J (a;) 
is due to the fact that the bath has to "resolve" the structure of the impurity, formally again the 
factor sin 2 (k • L). Furthermore, as already pointed out, no "ad hoc" cutoff frequency u c needs to 
be inserted but appears naturally in the decaying exponential of variance a in (129]) . 

Figj2] shows clearly that the impurity maintains much of its coherence at long times. Such 
survival is due to the above mentioned suppressed effect of soft modes, which are responsible of the 
long time behavior of ro(i), and is more pronounced when the environment consists of a condensate 
than in the case of a reservoir consisting of free bosons. This can be intuitively described in terms 
of greater "stiffness" of the condensate whose Bogoliubov modes are less displaced by the coupled 
impurity. The condensate is even able to give some coherence back to the impurity, since T^{t) is 
not monotonic in time. Oscillations of coherence in spin-boson systems were predicted in [26] (and 
even earlier, in a different context, in |35j). 

We can distinguish three stages in the dynamics of the r 's. In the first stage r (t) oc t 2 , as 
can be easily seen from a series expansion of (1291) . This very short stage, shown in the inset of 
FigfSJ corresponds to coherent dynamics. The second stage corresponds to a Markovian behavior, 
i.e., r (t) oc t, and lasts a few tens of microseconds. Finally, in the third stage T (t) saturates to a 
stationary value. This behaviour calls for particular caution in treating an environment of (free or 




(30) 



k 
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Figure 3. T\(t) (dashed line), Tzif) (dotted line), and 2To{t) (solid line) vs. time for a pair 
of impurity atoms at a distance 2D — 4L (see text), immersed in a condensate (left) and in an 
environment of free bosons (right) in three dimensions. 



interacting) bosons as a Markovian reservoir for atomic impurities immersed in it, which is clearly 
not the case in the present situation. 



4-2. Collective decoherence of two impurities 

Decoherence of quantum systems in a common environment is characterized by collective 
decoherence. It is well known that two spins interacting with the same bosonic reservoir with 
a spin - boson interaction Hamiltonian like the one discussed in this paper show sub - and super - 
decoherence [25]. In simple words the decoherence rate of the two spins is not simply 2Fo(t) but, 
according to the initial state of the spins, much smaller or larger. In this final section of the present 
manuscript we analyse the specific features of collective decoherence in our system. 

For two pseudospins, three decoherence parameters appear in the density matrix elements 
independently of the exact form of the impurities' state. One is T (t) and appears in elements such 
as po,0;0,i(£), Po,i-,i,i{t) , etc., which corresponds to individual dephasing of each impurity atom; two 
more parameters Ti{t) and T 2 (t) appear in elements such as |po,0;i,i (*) I = exp [— Ti(t)] |po,0;i,i(0)| 
and |po,i;i,o(£)| — ex P [~ T 2 (t)] |po,i;i,o(0)|, and corresponds to decay of the coherences between states 
with the particles in the same or in the opposite side, respectively, of the double well. For two 
pseudospins at distance 2D = 4L, these two parameters are 



ri(t) = r {0 , 0M1>1} (t) = ^ 2 ^ 



i 



cos iH 



El 



coth 



m k + q 2 



n3U 



32^ B n X)(l' u kl ~ Kl) 2 e 



-k 2 <T 2 /2 



^n 2 ft 
El 



sin 2 (k ■ L) cos 2 (k • D) 



(31) 



r 2 (t) = r {1 , 0M0>1} (t) = ft 2 £ 



(1 



cos iN 



El 



coth 



O 1 - O 1 



L,k 
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Figure 4. 5 c (t) (dashed line) and S n ' l '(t) (solid line) vs. time for a pair of impurity atoms in a 
three-dimensional environment. The inset shows S(t) for very short times < t < 2 /is. 



32g 2 AB n (Kl - Kl) 2 e 



-k 2 a 2 /2 



sin 



2 E± 
2h ' 



El 



sin 2 (k • L) sin 2 (k • D) 



Calculations similar to those performed for Tq give for a condensate environment 



(32) 



2g 2 AB n o 



dk k\^° 2 ^ ; in2 ^ , 

E k (e k + 2g B n ) 



sin2kL sm2kD sm2k(L + D) sin 2 HZ) — L) 
x I 2 - 2 — , h 2- 




2k{L + D) 



2k[D - L) 



r c 2 (t) = /I* JbV*^ 

2V ' 7T 2 JO 



E^ (e k + 25f B n ) 



sin2A;L sin2kD sm2k(L + D) sin2k(D - L)' 
x 12-2—^ 2 - 1 — - + 



2r (t) + 5 c (t) 



2k(L + D) 



2k(D - L) 



(33) 



(34) 



In the above equations it is easy to identify the term S c (t) which quantifies the deviation to the 
dechoherence exponent 2Fo typical of the decoherence of two impurities interacting with independent 
environments. Note that while To depends only on L i.e. on the spatial size of the double well, S 
depends non trivially on L ± D i.e. on the distance between the impurities of different wells. As 
before the special case of a bath of noninteracting bosons r^*'(t), r 2 '*'(£) are obtained from the 
above equations ( 1331) simply imposing g B = and E^ = e^. 

As in the case of single impurity decoherence the impurities do not loose all their coherence: 
Ti and T 2 saturate to a stationary value that can be varied with the help of Feshbach resonances. 
Furthermore Figj3] shows that in a system of two impurities coherence oscillations appear, both 
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Figure 5. Ti(t) (top) and ^(i) (bottom) vs. time for a pair of impurity atoms interacting with 
a bosonic condensate (left) and with free bosons (right) in three dimensions for different distances 
between the impurities: 2D — 8L (dash-dotted line), 2D = 16L (solid line), and 2D = AOL (dashed 
line); 2To(t) (dotted line) is also shown for comparison. 



for interacting and non-interacting bosons in the environment (even more pronounced oscillation 
are shown in Fig. [5]). Such coherence revival is due to the collective nature of the coupling, as 
quantified by S c (t) (5 n '*'(t) for free bosons). As shown in FigJUalso the 5(i)'s are characterized by 
three different time scales comparable to those analysed for r (t). In the first stage the difference 
\S(t)\ is negligible, since the presence of each impurity cannot have modified yet the environment 
seen by the other one; in the second stage, corresponding to the Markovian dynamics, the difference 
\S(t)\ steadily grows up; and in the third stage it decreases, reaching a stationary value. 

For a pair of impurities we observe super - and sub - decoherence, however with a peculiarity 
which is characteristic of the system here considered. Indeed we observe sub - decoherence in 
Ti = r{o,o},{i,i} and super - decoherence with T 2 = r{i >0 },{o,i}9 at variance with what one observes 
in a standard spin boson model, where their role would be exchanged [25]. This different behaviour 
is due to the particular configuration of our system: Ti gets contribution from superpositions of 
the states 1 , 0) and |1, 1), where the atoms sit in wells with identical distance, whereas the states 
|0, 1) and 1 1 , 0) , contributing to T 2 , correspond to atoms sitting in wells with different separations. 

Further insight on the features of the collective decoherence is gained by considering the 
decoherence of impurities sitting in sites which are at a larger distance than 2D = 4L = 600 nm. 
In Figj5] we plot the decoherence exponents for impurities trapped in lattice sites at distances 
2D = 8L, 2D = 16L, and 2D = 40L respectively. These plots suggest the following picture: initially 
the impurities decohere independently, as if they were each immersed in its own environment; at 
some later time, the environment correlations due to the impurities act back on them and give rise 
to oscillating deviations from 2r (t). The onset time of these oscillations depends on the separation: 
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the larger the separation, the later the onset. On the other hand, the correlations become weaker 
as the distance increases and the oscillations become consequently smaller in amplitude. At large 
separation (here, approximately 40L), the parameters I\ and T2 are hardly discernible from 2Fo, 
since the environment correlations induced by the impurities vanish. Similar features in a related 
context are reported in [36]. In summary Ti(t) and T^t) also prompt information about the 
correlation length of the environment. 



4-3. Decoherence in one dimension 



Finally, we examine the decoherence process in a one-dimensional condensate. Since, as previously 
discussed, the spectral density ( 1301) is superohmic for an interacting gas, but subohmic for a free Bose 
gas, we expect qualitative different results for the two cases, in contrast to the three-dimensional 
case. The decay exponents in one dimension ^(t) become 



lo(t) = 
for one impurity and 



4g 2 AB n o 



7T 



dk 

00 



sin^ ft 



E k (e k + 2g B n ) 



sin kL 



(35) 



4g 2 AB n o 



71 



dk 



2 7o (t) - 6 c (t) 



71 



dk 



= 2 lo (t) + S c (t) 



-fc 2 a 2 /2_ 



sin z f£i 



E k (e k + 2g B n ) 



E k (e k + 2g B n ) 



sin (kL) cos 2 (kD) 



sin (kL) sin 2 (kD) 



(36) 



(37) 



for two impurities in a condensate. The behaviour of these parameters critically depends on the 
nature of the environment, see FiglHl In particular, decoherence in a one-dimensional sample of free 
bosons results Markovian, in agreement with the naive expectation, due to its subohmic spectral 
density. 



5. Conclusions 



We have shown how a system of impurity atoms trapped in an array of double wells, interacting with 
a cold atomic gas, is described, in a suitable regime, by a spin - boson hamiltonian. The specific 
nature of our system, in which the pseudospins, associated with the presence of an impurity in the 
right/left well of each site, have a spatial dimension introduces peculiar features in the decoherence 
of a single impurity as well as in the collective decoherence, with the persistence of coherence at long 
times, the presence of coherence oscillations and counterintuitive super / sub decoherent states. 

We have shown in particular that for a three-dimensional bath one never has a Markovian 
behaviour. A one-dimensional bath is in this respect more interesting since one can go from a 
non-Markovian to a Markovian behaviour just by tuning the interaction of the bath. 
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Figure 6. 7i(i) (dashed line), 72(f) (dotted line), and 270(f) (solid line) vs. time for a pair of 
impurity atoms immersed in a condensate (left) and in an environment of free bosons (right) in one 
dimension. The separation between two impurity atoms is 2D = AL. 

As a final comment we would like to say a few words about the role of the quadratic terms in 
the Bogoliubov operators which we have neglected in our derivation of Hamiltonian (114j) . Although 
a detailed study of their effects is beyond the scope of the present article, we would like to point 
out that their effects are negligible with respect to the linear terms we have analyzed in the present 
manuscript. One can show that their inclusion amounts to taking into account elastic scattering of 
Bogoliubov particles, which is simply responsible of an energy shift, inelastic scattering processes 
and Bogoliubov pair creation and annihilation. In these two latter additional terms the length of 
wave vectors k that can play some role in the impurities' dynamics is limited from below by the 
finite size of the condensate and from above by cutoff parameter er -1 . It can be shown that, in this 
frequency range, the coupling constants of the neglected processes are, for the values of parameters 
assumed in our analysis, three orders of magnitude smaller than the coupling constants hQ l n]l of 
the linear terms. As a consequence, a rough estimate leads us to suppose that any possible relevant 
effect of the quadratic terms in the Hamiltonian would become apparent at time scales that are 
three orders of magnitude larger than those examined in this article. 
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Appendix A. Disentangling the time evolution operator 



The factorization of the time evolution operator U(t) = exp —iHt/H is often an impossible task. 
When the Hamiltonian contains operators forming a Lie algebra the transformation of U(t) into a 
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product of simpler exponential operators is however possible in some cases [37]. Here we show a 
practical way to transform U(t), which we write as 



U(t) 



ex p -j-J2 E ^At ex P J2 ^E A k(*)^ + « k (*)l 4 



x exp 



u R (t) 



(A.l) 



where £7r(£) is to be determined, as well as the quantities A l k (t), B^(t), a k (£) and flk(t). Since at t = 
the time evolution operator {7 reduces to the identity operator, A k (0) = £> k (0) = Ac(0) = « k (0) = 0. 
All unknown quantities can be found with the help of the relation 



H = ih 



dU{t)/dt U~\t) 



which holds for any time-independent Hamiltonian and of the relation 



e x Ye 



-x 



Y + 



X,Y 



+ 



X, 



X,Y 



+ 



1 r 



X, 



X, 



X,Y 



+ ... 



(A.2) 



(A.3) 



for arbitrary operators X and Y. After inserting the expression (1 A. II) for the time evolution operator 
U(t) in the right-hand-side of flA.21) . a comparison with the Hamiltonian ffl8l) leads to the expressions 



a k (t) 



a i n R k ii L k 



2£ k 



(l - e^) , = a k (t) 



for -£>(£), and Pit), and to the differential equation 



(A.4) 
(A.5) 



d_ 

dt' 



-Mt)= -e(e^(^ + 4w) \^Amai + a k (t)ju R (t) 

for the unknown exponential operator Unit), which we write as 



U R (t) = exp 



(A.6) 



(A.7) 



A comparison with (IA.60 gives 

ri{t) = Bl(t)A{(t), e k (t) = M)a k (t), Hit) = B£(f)«k(*) + &(*K(f) (A.8) 



i.e., 
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— %h 



Eij (nk,k + ^L, k ) Kk + n£ 



-iE k t/h 



ill 



4£ k 

- «U) E Kk + 4: 



/ + — (l 



E, 



ih 



Lib / 



-iE k t/H 



(A.9) 
(A.10) 
(A.ll) 



Moreover, using Glauber's relation 



exp 



E#k4 

. k 


exp 


- E f kC k 
k 


= exp 


I] (#k4 - ^k4) 

. k 


exp 


~E \9k\ 
- k 



(A.12) 



the two exponentials linear in Bogoliubov operators can be merged into 

exp E(E^)^ + «k(0)4 exp -E(E 5 k(^ + /5k(t))c k 

= exp | E (e 4(*)^i + «k(*)) 4 - E (e 4*(*)^ + 4 j 
E fe + ok(t)l IE 4*(*)<^ + 



(A.13) 



x exp ^ - 



(A. 14) 



and the contribution of the last exponential can be included in Un(t). Performing some 
commutations where it is possible, the time evolution operator becomes 



U(t) = exp 



exp 



k \ ij i 



x exp 



x exp 



r E EktfAJ 
1 k 

E (e Aim + a k (*)l 4 - E (E + ) c : 



(A.15) 



Finally, the exponential operators that do not contain bath operators commute, so the time evolution 
operator can be further modified into the final form (jl~9]) . 

Appendix B. Derivation of the dynamics of the impurities 

The action of £/ k) z>(t) on a pure state of the whole system is 
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UKD(t)\{n p })({m p }\ ®p k Ul D (t) = \{n p })({m p }\ ® 



cxp 



exp 



(B.l) 



and the density matrix elements P{n p },{m p }(t) of the impurities are obtained by tracing over the 
bath, 

P{n P },{m P }(t) = exp {iQ{n p },{ mp }(t)}exp {^{n p },{m p }(t)} p {riph{mp} (0) 

x (K}| l[Tr B>lc {u k , D (t)\U P }){{m p }\ ® Pk Ul >D (t)} \{m p }) (B.2) 

k 

where Tr^ k denotes the trace over each Bogoliubov mode of the environment and the phases 



@{n p },{m p }(t) = h 2 J2 



k 4££ 



x 



E » Kk - «U) Kk - [(-i) n ^ - (-i)^] (b.3) 



= ^ 2 E ^#^E Kk + n&c) E Kk - «U) (m - ^) (b.4) 



k ^ 



come from the unitary operators in ([191) . Performing ciclic permutation inside the trace and using 
the identity exp(M) exp(iV) = exp(M + N) exp[M, N]/2, which holds for operators M and N that 
commute with their commutator, the trace Tr^k in (1251) becomes 



exp 



x Tr 



B.kjexp 2X)(n i -m i )(4(t)4-^(t)c k j /> k 

= exp {tA {nph{mp} (t)} Tr Bjk jexp 2 ^ (n, - m 4 ) (^(^4 ~ 4* 
The trace over the thermal bath of the displacement operators is well known [25], 



Pk (B.5) 



Tr 



B,k 



exp {# k c k - # k c k } p k 



exp 



l2k| 2 f ,^k 
coth 



(B.6) 



where /3 = (ksT) 1 , and leads to equation (125]) . 
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Appendix C. The coupling constant in a deep optical lattice 

In a deep optical lattice the ground state wave functions of each well can be approximated with 
those of harmonic oscillators, 



exp 



(x - x i:N ) 2 (y - y iiN ) 2 (z - z i>N f 



2Xn 



22/o 2 



Izl 



(c.i; 



Here N = L,R, and xq = ^Jh/(mu x ), y = ^Jh/(mujy), and z = ^Jh/(mu z ), where the u's are the 
trapping frequencies of the harmonic trap approximating the lattice potential at bottom of L and 
R wells of the lattice site i. The coupling frequencies (fTTj) of the spin-boson model then become 



n 



_ 9ABy/n 

- k ~~ n 
h 



(|«k| - \v k \) d x \(p i)L (-x)\ e 



(|«k| - Kl)e 



2 ikx 



-fc 2 cr 2 /4 ik x Xi^ n 



n = L,R 



(C.2) 



having assumed identical confinement in the three directions, a = x = y = z . 
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